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The modal interpretation of quantum mechanics allows one to keep the standard classical 
definition of realism intact. That is, variables have a definite status for all time and a measurement 
only tells us which value it had. However, at present modal dynamics are only applicable to 
situations that are described in the orthodox theory by projective measures. In this paper we 
extend modal dynamics to include positive operator measures (POMs). That is, for example, 
rather than using a complete set of orthogonal projectors, we can use an overcomplete set of 
nonorthogonal projectors. We derive the conditions under which Bell's stochastic modal dynamics 
for projective measures reduce to deterministic dynamics, showing (incidentally) that Brown and 
Hiley's generalization of Bohmian mechanics [quant-ph/0005026, (2000)] cannot be thus derived. 
We then show how deterministic dynamics for positive operators can also be derived. As a simple 
case, we consider a Harmonic oscillator, and the overcomplete set of coherent state projectors 
(i.e. the Husimi POM). We show that the modal dynamics for this POM in the classical limit 
correspond to the classical dynamics, even for the nonclassical number state |n). This is in contrast 
to the Bohmian dynamics, which for energy eigenstates, the dynamics are always non-classical. . 



I. INTRODUCTION 

It has been almost eighty years since quantum theory 
emerged as a complete theory, yet fundamental debates 
still occur over its interpretation. These debates usu- 
ally center around what has been called the measurement 
problem, which we argue is a twofold problem. 

In quantum mechanics a complex vector which 
belongs to a Hilbert space Tt , is used to describe the state 
of a system, and its evolution is given by the Schrodingcr 
equation. In the orthodox interpretation it is postulated 
that a system only has a definite value for an observ- 
able A (energy, position etc.) if the state \^(t)) is an 
eigenstate, \^ n (t)), of A, in which case the value of A 
is the eigenvalue a n associated with \^ n (t)). However, 
since the Schrodingcr equations in linear it is possible 
for the system to be in a state which is a linear super- 
position of these eigenstates: \^f(t)) = S n Cn(i)|$ n (i)). 
That is, the system can have two or more values of an 
observable at once or, we say, the value is not defined. 
To explain why only one of the values is obtained in a 
measurement we need to introduced an extra dynamical 
equation, the reduction equation. This equation for a 
measurement of observable A has the effect of collapsing 
the state, — > \^f n (t)}, thereby defining the value of 
the observable. Thus the standard view of reality (that 
variables have well defined values even when they are not 
observed) is lost. This is what we call the first problem 
of the measurement problem under the orthodox inter- 
pretation. 

The second problem is concerned with choices: what 
observable is measured, and where the reduction occurs 
(the Heisenberg cut [1]). In the orthodox interpretation 



there is a classical world, which we live in, existing out- 
side of the quantum world and allows us to define an 
apparatus which chooses the observable to be measured 
(the set of eigenstates to be collapsed into). The problem 
is, at what point do we place the distinction between sys- 
tem and apparatus? An apparatus is made of particles 
just like the system. Perhaps it is only an intelligent ob- 
server that collapses the wavefunction. But it is not clear 
how intelligence or consciousness can influence physics. 
Alternatively it can be argued that quantum mechanics 
should extend up the von Neumann chain [2] to include 
the entire universe, so \^(t)) now labels the state of the 
universe. If this is the case then what is external to the 
universe which makes the quantum measurements? 

In this paper we consider the modal interpretation of 
quantum mechanics [3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 
15]. This interpretation's central goal is to keep the stan- 
dard definition of reality intact, that observables have 
values even when not observed. In this interpretation to 
explain why only one value of A is actual for all time 
we introduce an extra state, the property state, \^ n (t)}. 
This state evolves stochastically with time (jumping be- 
tween different values of n) and selects the present value 
of A, a n , from the set of possible values {a n }. While 
\^(t)) in effect chooses the weights for this stochastic 
jumping; that is, it acts as a guiding state. We define 
the unnormalizcd (signified by a tilde) property state as 

!*»(*)> = *»(*)!*(*)>, (LI) 

where 7t n (t) is a projector (acting in the total Hilbert 
space of the universe 7i U ni) and satisfies ^2 n ^ n (t) = 1. 
The set of projectors, {ir n (t)} are labelled the preferred 
projective measure. The normalized property state is 
then defined as 

|*n(*)}=*«(*)l*(*)>/^ (1-2) 
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where A is a normalization constant. Note the set of 
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property states depends on both the preferred projec- 
tive measure and the guiding state. The guiding state 
\^(t)) (as in the orthodox interpretation) is found from 
the Schrodingcr equation 

dt|*(t)) = -iH uni (t)\*(t)}, (1.3) 

where H un [(t) is the Hamiltonian of the universe. 

The reason why only one value for A is obtained upon 
a measurement of A at time t, is we choose the property 
states (i.e. choose fr n (t)) so that A can be written as 
J2 n a n^n(t)- Now if |^„(<)) is the property state at time 
t, then A has the value a n at t and if A is measured, this 
value is revealed. The probability that universe will be 
in the n th property state at time t is given by 

P n {t) = (*(f)|* B (t)|*(f)) = (tt n (f)|tt n (t)>, (1.4) 

the standard Born probability. Thus the average of A 
over all possible property states will agree with (A) found 
by the orthodox theory. However, unlike the orthodox in- 
terpretation this probability does not refer to the proba- 
bility of observing results a„ at time t but to the proba- 
bility of A having value a n at time t. 

For a given set of property states, there will be more 
than one observable, which can be assigned a definite 
value. Such and observable should, and from now on 
will, be referred to as a property or beable (after Bell 
[3]). Our notation for a property is 

Z = {(z n ,n n (t))} (1.5) 

where z n is the value of the property, which could be a 
real number, or a complex number, a statement (yes/no) 
or even a string of numbers. When the set {z n } are real 
(complex) numbers the property can be represented by a 
Hermitian (normal) operator 

Z{t)=J2 Z nKn(t). (1.6) 

n 

Thus, in this interpretation, the first problem of the 
measurement problem is solved, at least for measure- 
ments of the preferred observable (property), as a mea- 
surement only tells us which value was possessed. For 
measurements of other observables see the discussion by 
Bohm for their interpretation [16]. However, the prob- 
lem of choice still remains as we can choose a different 
preferred projective measure, {J n (f)}, which will give a 
different group of observables property status. This prob- 
lem for modal interpretations has not been resolved in a 
wholly satisfactory way, and may never be. 

The problem of choice is illustrated by the many vari- 
ants of the modal interpretation. In the beable vari- 
ant [3, 4, 5, 6, 7] choice of the preferred projective 
measure is fundamental. That is, it is made inde- 
pendently of the guiding state. The problem is that 
many choices are viable [6]. In other modal theories 
[8, 9, 10, 11, 12] the preferred projective measure depends 



on |^(t)). These have been labelled by Bacciagaluppi 
and Dickson [9] as the atomic version of the modal inter- 
pretation. One assumes a preferred factorization of the 
universe TL a (g) Tip (g) . . . ® Ti^ . (This seems to be necessary 
to avoid a Kochen-Specker type contradiction [17]). In 
this preferred factorization the spectral resolution of each 
subsystems p" ed (t) = ^ n . w ni tc%. (i) defines the preferred 
projective measure for the universe as 

{U*) = ^(*)<^)--<(*)}, (1-7) 

where n = (m ...,n z ). This may seem like it has rec- 
tified the problem of choice but it still remains here in 
the choice of the preferred factorization. Another prob- 
lem of choice can emerge in the atomic variant this being 
that for some guiding states there are many viable algo- 
rithms for defining the preferred projector. Note other 
variants of the modal interpretation [6, 13, 15] which try 
to answer choice also encounter the above problems. 

In this paper we do not propose to solve the problem 
of choice. To the contrary, we show that there is actu- 
ally more choices to be made than what was previously 
realized. Specifically, it is not necessary to restrict the 
property to be of the form of Eq. (1.5), based upon or- 
thogonal projectors. The extra choice is motivated by the 
fact that in the orthodox interpretation of quantum me- 
chanics we can extend the theory of measurement from 
that using projective measures to that using positive op- 
erator measures (POMs). The former measurement (of 
observable Z) is described by a set of pairs as in Eq. (1.5), 
but here the z' n s would be interpreted as measurement 
results. The latter measurement (of observable Z) is de- 
scribed by a set of pairs 

Z = {(z n ,F n (t))}, (1.8) 

where F n (t) is referred to as a POM element [18] or ef- 
fect [19]. These effects are positive and complete, with 
^2 n F n (t) = 1. In this paper, we show that we can de- 
velop modal dynamics for preferred observables of the 
form in Eq. (1.8). 

To agree with orthodox theory, our modal dynamics 
for POMs should generate the usual probability formula 
for POMs, which is 

P n (t) = (*(f)|F n (f)|tf(t)). (1.9) 

Note that in general the number of possible results 
(Y]„ 1) can be greater than the dimension of 7i U ni- Some 
examples of POM measurements include informationally 
complete POMs [20, 21] and the Husimi or Q-function 
POM [22, 23]. In the orthodox theory the reduction 
equation for POM- type measurement is — > M n |\P) 
where M n is referred to as the measurement operator 
and is defined such that F n {t) = M\M n [19]. 

Since the modal dynamics presented in Refs. [3, 4, 
5, 7, 8, 9] assume that the property states are distin- 
guishable, non-orthogonal states are not allowed in their 
analysis. In this paper we use Naimark's theorem [18] to 
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extend the dynamics to include POMs. We begin in Sec. 
II by extending the stochastic dynamics introduced by 
Bell [3]. In Sec. Ill we consider the limits in which deter- 
ministic dynamics can be derived from the stochastic dy- 
namics, for both the projective and the POM case. Along 
the way we show that Brown and Hilcy's [25] generaliza- 
tion of Bohmian mechanics cannot be derived from Bell's 
modal dynamics. In Sec. IV we illustrate the POM case 
with an example of a universe consisting of only a sin- 
gle Harmonic oscillator. We show that when the Husimi 
POM (the overcomplete set of coherent state projectors), 
is used, the modal dynamics is consistent with classical 
mechanics in the classical (large excitation) limit. How- 
ever when the modal dynamics is used to describe posi- 
tion property states (i.e. Bohmian mechanics) there is no 
recovery of classical behaviour in the classical limit. This 
illustrates the potential usefulness of modal dynamics of 
POMs. 



II. DYNAMICS FOR MODAL 
INTERPRETATIONS 

In this section we first briefly outline the modal dy- 
namics for the projective case, which is basically a re- 
production of the results presented in Refs. [3, 4, 5] and 
generalized in Refs. [7, 8, 9] to include time-dependent 
projectors. Secondly we use Naimark's theorem [18] to 
extend the dynamics to include positive operators. 



A. The Projective Case 

Before reproducing the standard results for the projec- 
tor case we outline the general method used to describe a 
stochastic process, that jumps between N distinct states. 
We define P n (t) as the probability that the system is in 
the n th state at time t. Assuming a Markovian process, 
by which we mean that the probability of being in state 
m at time t + dt only depends on the state at time t. we 
can write the master equation as 

d t P n (t) = Y^[ T nm(t)P n {t) ~ T mn P n {t)], (2.1) 
m 

where T nrn is defined as 

T nm (t) = [P nm (t + dt, t) - P nm (t, t)]/dt. (2.2) 

Here P n m(t + dt,t) is a conditional probability and is 
read as the probability of the system being in state n at 
time t + dt given it was in state m at time t. From this 



definition it follows that ^ n T, 



n 1 nm 



0. 



For distinct states the conditional probability P nm (t, t) 
must be for n ^ m. This allows us to interpret T nm (t)dt 
for n in as the transition probability, so we call T nm (t) 
the transition rates. For n = m, P nn (t, t) = 1 and T nn 
(which is negative) is a measure of the rate at which state 
n losses probability. 



Defining a probability current J nm (t) as 

results in J nm (t) = —J m n(t) and allows us to rewrite the 
probability master equation as 



dtPn(t) =J2 J nm{t). 



(2.4) 



Given J nm {t) and P n (t), there are many possible transi- 
tion rates satisfying Eq. (2.4). One solution, chosen by 
Bell [3] is as follows. 
For J nm (t) < 0, 



^mn (^) ~~ Jnm (fy j Pn (^) 3 

and for J nm (t) > 

T nrn (t) — J nm [tjfP m (t ) , 

T mn (t) = 0. 



(2.5) 
(2.6) 

(2.7) 
(2.8) 



This is only one of the infinitely many solutions. These 
are found by adding an extra term, T® m (t), to T nm (t), 
where T® m (t) is constrained only by 

T° m (t)P m (t) - l£ n (t)P„(f) = 0. (2.9) 

To make the link with quantum mechanics we say that 
the distinct states are the property states {|<I' n (£))}, 
and the set of properties which have definite values form 
the group Q, defined by elements 



Z = ]Pz„7r„(i), 



(2.10) 



with the value of property Z being the corresponding z n . 
The evolution of these values (jumping between z n ) is de- 
termined by the rates T nm (t), which themselves depend 
upon J nm {t). 

By taking the time derivative of Eq. (1.4) we obtain 
the differential equation 

d t P n (t) = 2lm[(^(t)\n n (t)H uni (t)\^(t))} 

+(*(t)|dt[7r n (t)]|*(t)), (2.11) 

where we have used Eq. (1.3), the Schrodinger equation. 
This can be simplified by defining the Hcrmitian opera- 
tor, R(t), for which 



d t K(t) = -i[R{t),Tt n {t)]. 
This allows us to rewrite Eq. (2.11) as 



(2.12) 



d t P n (t) = 2lm{(^(t)\Tr n {t)[H uni {t) - A(t)]|tf (*))}. 

(2.13) 

Comparing this with Eq. (2.4) and using the fact that 
X] m 7rm(t) = 1, one possible probability current (that 
was chosen by Bell [3]) is 

J nm (t) = 2lm{(^(t)\7r n {t)[H uni (t) - R(t)]it m (t)\9(t)}}. 

(2.14) 
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Note that this is only one of infinitely many possible cur- 
rents, as we can add any current J^it) to J nm {t) which 
satisfies J2m Jnm = 0> to Si ye a valid probability current. 
For the purposes of this paper we only consider the sim- 
ple solutions (not containing the extra T® m (t) and J r ® n (t) 
terms). For a discussion about these solution see [5] and 
[8]. 

B. The Positive Operator Measure Case 

A typical example of a POM is the observable 

Z = {(z n ,^\z n }(z n \)}, (2.15) 

where {|z n )} forms an ovcrcompletc basis in 7i U ni, with 
\z n ) (z n \ — IB. If the universe was in \z n ), then a 
measurement of this observable could with finite proba- 
bility yield any of the values {z n }. This clearly disagrees 
with the above description of modal dynamics. Thus we 
cannot treat non-orthogonal states in the same manner 
as orthogonal states. 

POMs are not always generated from non-orthogonal 
states. Given any set of projectors {II„(t)} in a larger 
Hilbert space K, = 7i un i <8> Haux where 7i aU x is some aux- 
iliary Hilbert space, it is well know that a set of POM 
elements can always be found by [18] 

F n {t) = Tr )n„(*)] (2.16) 

where p aux is a state in 7i a ux- For simplicity we define it 
as paux = What is perhaps less well known is that 

given a POM {F n (t)} it is always possible to generate 
a projective measure {II n (t)} in a larger Hilbert space 
tl un i <g> H a ux, where dim(7i aux ) is not necessarily equal 
to X) n l. This is called Naimark's theorem [18, 24]. It 
basically says that we can define a projector TL n (t) such 
that 

Tr Ha JMt)){*(t)\F n (t)] 

= H:jc[|*(t)><*(t)| ® |0>(0|n n (i)], (2.17) 

for all |*(t)) G 7Y U ni and for n = 1, N where N is the 
number of POM elements in the POM. To work out the 
set {II„(t)} it is necessary to introduce another projector 
tlN+i(t) 7 such that 

JV+l 

uni+aux; 

(2.18) 

n 

and 

n„(t)ft m (t) = fl n (t)6 nm , (2.19) 

is satisfied for n, m = 1, N + 1. The set of projectors 
in this enlarged Hilbert space is called the Naimark ex- 
tension of F n (t) [18]. A worked example of this is shown 
later in this section. 



We now propose that to calculate modal dynamics for 
POMs, the preferred observable (determined by the pre- 
ferred POM) defined by Eq. (1.8) becomes a property in 
/C. The set of properties which can have definite values 
is now defined by 

Z = {(z n ,fl n (t))}, (2.20) 

and if z n is a number this is represented by the operator 
Z = J2n +1 z rS^ n {t). The value zjv+i is arbitrary. Note 
we expect that the preferred projective measure {n ra (i)} 
may define more then one POM (set of POM elements) 
as preferred. The guiding state becomes 

|*(f)> = |*(*)>®l$> (2.21) 

where is still the solution to the Schrodingcr 

equation (1.3). With this guiding state we can rewrite 
Eq. (1.9) as 

P n (t) = ($(t)|H„(t)|$(i)), (2.22) 

which forces Pjv+i(t) = for all time as the projector 
njv+i(t) projects into the null space of | $(£))• The prop- 
erty states are defined in K, as 

]$„(*)) =II„(t)|$( i ))/ViV 7 , (2.23) 

which now form an orthogonal (distinguishable) set and 
N' is the new normalization constant. Thus the stan- 
dard analysis of modal dynamics now applies (Sec. II A). 
However, the property state in this extended Hilbert 
space is, in general, an entangled state (between the uni- 
verse and the auxiliary system). This may raise inter- 
pretational difficulties, but we will leave this question for 
later work and just treat the above as a purely mathe- 
matical procedure to deal with POMs. 

To work out the modal dynamics for POMs, as in the 
projector case, we need to find the probability current. 
Using Eq. (2.22) and Eq. (1.3), J nm (t) is 

J nm (t) = 2Im{($(i)|n„(t)[^ uni (i)(8)i aux -E / (t)] 

xn m (i)|$(i)>}, (2.24) 

where R'(t) is Hcrmitian and defines the evolution of the 
projectors tl n (t), by 

d t fl n = -i[R'{t),il n }. (2.25) 

To illustrate the above we consider a simple exam- 
ple, a universe consisting of only a spin- 1/2 system [18], 
and consider the observable Z defined by the POM el- 
ements F n = £\z n )(z n \, with n being 1, 2, and 3 and 
z n = exp(i27m/3). The states \z n ) are defined by 




(2.26) 



where |2) and |1) are the eigenstates of a z (the Pauli spin 
matrix). In the Bloch sphere these states all lay in the 
x — y plane with an angular separation of 2ir/3. 
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Using Naimark's theorem we extend this 2-dimcnsional 
Hilbert space to a 4-dimensional Hilbert space, where it 

I 

til 

n 2 
n 3 
n 4 



The two states \(j>) and \<f>') form a set of orthogonal basis 
states of the auxiliary Hilbert space and II4 is an added 
projector needed to complete the set of projectors. 

The modal dynamics for these states is formulated as 
follows. By Eq. (2.21) and Eq. (2.23) the three possi- 
ble property states are {| = n n |0i) |*(t))/v / iV 7 }, 
where n = 1,2,3. The fourth projector is not included 
as for all possible states it projects into the null space 
of |$(t)). The observable Z = {(z n ,F n )} becomes the 
property Z = {(z n ,tl n )} (or Z = X) re z„n„ since for all 
fly Zfi IS cl number) , and for this example the possible val- 
ues are; z\ — e ?27r / 3 , z 2 = e ?47r / 3 , Z3 = 1. The stochastic 
evolution between these values is then found from the 
transition rates T nm (t), which in turn require a specifi- 
cation of J nm (t). For simplicity, and to illustrate the sig- 
nificance of enlarging the Hilbert space, we assume that 
there is no evolution, that is J nm (t) = T nm (t) = 0. Thus 
once we first assign a property state (based on the initial 
probability distribution) it remains in this property state 
for all time. Now we assume that the initial state of the 
universe is in one of the three non-orthogonal states, say 
|*(0)) = \zi). Unlike in the projector (orthogonal) case 
we expect there to be some probability for the universe to 
be in any of the 3 possible property states | $„(£)). This 
does occur as in general n„|0)|zi) is non-zero for all n. 



III. CONTINUOUS MODAL DYNAMICS 

In this section we investigate the continuum limit of 
modal dynamics. This has been previously done by Sud- 
bcry and Vink in Rcfs. [4] and [5] respectively, where 
it was shown that Bohmian mechanics can be obtained 
by choosing the appropriate preferred projective measure 
(property states). In this section we briefly outline their 
work, then extend it to show that the Brown and Hilcy 
generalization of Bohmian mechanics to include the mo- 
mentum representations [25] is not the continuum limit 
of Bell's modal dynamics. We also present an alterna- 
tive (we believe easier) method for finding the continuous 
trajectories, that works when the modal dynamics has a 
continuum limit. 



can be shown (using Eqs. (2.18) and (2.19)) that the four 
projectors are [18] 



(2.27) 

I 

A. Bohmian Mechanics 

In Bohmian mechanics [16] the preferred projective 
measure is the one associated with the position basis. 
That is, the property is the position of the system {Xj} 
(vector notation) and the wavefunction ty({xj}, t) — 
({xj}\^>(t)) is then interpreted as a field which guides 
the position in a non-classical way. In Bohm's original 
papers he showed that this non-classical behavior could 
be represent by an extra potential in the Hamiltonian- 
Jacobi equation, the quantum potential which depends 
on ¥({*,•},«). 

In this paper we do not introduce the quantum po- 
tential, but instead describe Bohm trajectories with ref- 
erence to a continuous probability current Jk({xj},t). 
We also consider [^(t)) to belong to a tensor product 
Tt x ® H', where TL X is the Hilbert space containing the 
position eigenstates Then the wavefunction be- 

comes a vector given by, 

m{x j },t)) = {{x j }\*(t))€H f . (3.1) 

This allows us to take into account a Hilbert space for the 
universe which is larger than that for the position of the 
system. Here we see that |{iEj})|V'({ :c j}) U) i s the contin- 
uous equivalent of our unnormalizcd property state, and 
the properties are the position operators {^0c}- With 
this property state we can define a probability density as 

Pi{x j },t) = $({x j },t)${{x j },t)) (3.2) 
which obeys the continuity equation 

d t P({ Xj },t) = -^d^Mix^t), (3.3) 
k 

where Jk({xj},t) depends on the form of H nn i(t). As 
in the modal case there is not a unique solution to 
Jk({xj},t) as for example in three dimensions we can 
add any vector field V x A({xj},t) to 3({xj},t) as 
V • V x A({xj},t) = 0. 

We define a velocity field Vk({xj},t) implicitly by 

M{xj},t) = v k ({x J },t)P({x J },t). (3.4) 

Bohmian trajectories for {xj}(t) are then defined by 

d t x k {t) = Vk({xj},t)\ Xk=Xh (ty (3.5) 



= A ® 10)^1 - V2/3|zi)(z 2 | ® \4>){<f>'\ - V2/3|z 2 )(zi| ® \<t>'){<t>\ + l/3|z 2 )(z 2 | <g> \tf>'){4>'\, 

= F 2 ® \4>}{(f>\ + V2/3|z 2 )(z 2 | ® \4>)(<j>'\ + V2/3|z 2 )(z 2 | ® \cf>')(cf>\ + l/3|z 2 )(z 2 | ® \<t>'}(4>'\, 

- F 3 ® \4>){cf>\ + y2/3|z 3 )(^ 2 | ® \4>}(<t>'\ + V2/3|z 2 )(z 3 | ® + 1/S\z 2 )(z 2 \ ® \<t>'}(4>'\, 

= i<»|^|-|z 2 )<z 2 |®|^|. 
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Probability enters only through the initial conditions, 
{xk(Q)}- This is only one of the possible continuous tra- 
jectories which satisfy Eq. (3.2). Other possibilities in- 
clude stochastic approaches; see Rcf. [5] and references 
within. It should be noted that so far nothing has speci- 
fied that {xj} must be positions, and in fact Brown and 
Hiley [25] develop a formalism where {xj} can be cither 
position or momentum. As examples, they consider a 
simple universe (a single 1-D particle) and derive dtx(t) 
and d t p(t) for a linear, quadratic and cubic potential. 

1. Bohmian mechanics as the continuous limit of discrete 
modal dynamics 

To demonstrate that the modal dynamics does give 
Bohmian mechanics as its continuum limit, first consider 



the Hamiltonian 

#uni(*) = A(t) + Y,B k (t)Y k +Y,Y k Bl{t) 

k k 

+£&(*)% 2 +£Wt(*)> ( 3 - 6 ) 

k k 

where A(i), B k (t) and Ck{t) arc arbitrary functions of 
the operators {Xk} and the rest of the universe, and {Y k } 
are the conjugate operators to {Xk}- That is, \Yj, Xk] = 
~iS jk - 

To calculate d t Xk(t) in Eq. (3.5) we need to calculate 
the velocity field Vk{{xj},t), which in turn requires cal- 
culation of Jk{{xj}, t). Taking the derivative of Eq. (3.2) 
and using the Schrodingcr equation (Eq. (1.3)) we can 
write 



d t P{{ Xj },t) = 2Im[(^(t)\{x j }}({x j }\H uni (t)Mt))} 
= -2^Re[(Vi({^,i)|{2Re[4(^ 

k 

-2£^Mte},i)|{2Re[^ 

k 

(3.7) 

which can be shown to be 

d t P({ Xj }, = d Xk {ite [$({*,-}, t)\{2Re[B k ({ Xj }, t)] - UMCk({xj}, t)]d Xh - 2id Xh t)]}$({xj}, t))] }. 

k 

(3.8) 

Comparing this with Eq. (3.3) gives 



Jk{{Xj},t) 



Re 



M{s 3 -},t)|{2Re[B fc ({x 3 -},t)]-^ (3.9) 

I 



For simplicity we consider the case when A(t) = V(X), where Rj(x,t) and Sj(x,t) are real functions and {\j}} 

B k (i) = 0, and C(t) = &(t) = 1/4M, describing for is an orthonormal basis set, which for example spans the 

example an electron in a 1-D potential, with the operator Hilbert space of the internal structure of the electron, 

nature of 1 and V signifying operation on the Hilbert Then it simplifies to 



space for the internal structure of the electron. For this 
example J(x,t) becomes 



1 

Jl 



d t x(t) 



EjR^x,t)d x [s 3 (x,t)} 



J(x,t) 
and thus 

dtx(t) = —Ira 



Im[(4>(x,t)\d x \i,(x,t))}, 



$(x,t)\d Xk $(x,t)) 



x—x(t) 



(iP(x,t)\iP{x,t)) 
To simplify this we can rewrite \ip(x,t)) as 

\i>(x,t)) =Y / Rj^,t)exp[iS j (x,t)]\j) 



(3.10) 



(3.11) 



(3.12) 



x—x{i) 



(3.13) 



To compare this to the modal dynamics defined in 
Sec. II we have to discretize X. In [5] and [4] this is 
done by defining a lattice of size and lattice sepa- 
ration e. Thus the values of the property X denoted 
x become x n = en, and the preferred projective mea- 
sure becomes {7r„ = |en)(en| (g> 1}. With this preferred 
projective measure the property state Eq. (1.2) becomes 
|*n(*)) = \eri)(en\iff(t)) = \en)\^ n (t)) where \ii n {t)) is an 
unnormalized state existing in TC . Using the results of 
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Sec. II A the probability current is 

J nm (t) = 2Im[(^ n {t)\(en\H uni \em)\^> m (t))}, (3.14) 
and the discretized version of the Hamiltonian is 
(en\H uni \em) 



— (<5n,m+l + ^n,m-l — 2<5„, m )/2Me 2 

+V(en)S ntm . (3.15) 



This gives 

J nm (^) 



— — : jlm[(-0n (t)\lpn-l(t))6 n . m+ i 

+(i> n (t)\i> n+1 (t))6 

n,m — 1J • 



(3.16) 



Taylor expanding \ip n+ i(t)) and |^ n _i(t)) gives 
1 



Me 



Im{(</>„(i)|A e [|V> n (*))]<S 



71,771+1 



■(^ n (t)\A e [\^ n (t))]6 n , m -i + 0(e)}, (3.17) 



where A e is the discretized version of a derivative. As in 
the continuous case we write \i/j n (t)} in terms of the real 
functions Sj(en,t) and Rj(en,t) which results in J„ m (i) 
becoming 

J nm (t) = — R 2 j (en) A e [Sj (en)]<5 n , m+ i 

i 

-^i? 3 2 (en)A e [S 3 -( £ n)],S nim _ 1 + 0(6)}. 
i 

(3.18) 

Since 1/e ^> 1, in the e — > limit (continuum limit) 
we can neglect the higher order terms in the above ex- 
pression for J nm (t). That is the only terms which will 
contribute are the transitions from m to m — 1 or m to 
m + 1. Because of this we can write 



J n{n+1) (t) =- -L^^^A.^M] (3.19) 

i 

= ^I>iMA £ [S,-(en)]. (3.20) 
If V i^A^en)] > (implies J n{n+1) {t) < and 



■k(n-i)(*) > 0) then by Eqs. (2.5) - (2.8), 

->-(n+l)nV') — 



Me£,J8(en) 



T( n -i) n (t) 



(3.21) 
(3.22) 



If "£.R^A e [Sj(en)] < (implies J n(n+1) {t) > and 
(*) < 0) then by Eqs. (2.5) - (2.8) 



T'(ri+l)n(^) 



(3.23) 



T, 



(n-l) 



»(*) 



ZjR^enjA^Sjien)} 



Me£,.ii?(en) 



These transition rates imply that in an interval dt the 
average displacement dx will be 



E[dx] 



e T( n+ i) n dt — eT(„_ 1 ) n df 



di + 0(e)rfi. (3.25) 



Provided Sj(en) and Rj(en) are continuous, the average 
-Efefa;] reduces to Eq. (3.13) as e — > 0. However, to show 
that the trajectories are smooth and deterministic from 
the initial {xk{0)} we also require that the dispersion 
E[dx 2 ] goes to zero in the continuum limit. This is the 
case since 



E[dx 2 } = e 2 T {n+1)n + e 2 dtT {n _ 1)n dt = 0(e)dt, (3.26) 



which goes to zero as e — > 0. 



2. Hamiltonians for which Bohmian mechanics is not the 
continuous limit of discrete modal dynamics 



The above demonstrates that in the continuum limit 
modal dynamics becomes Bohmian Mechanics. However, 
if we consider a Hamiltonian of the form 



H ur .; — kY 



V(X), 



(3.27) 



we find that the continuous limit of discrete (Bell-type) 
modal dynamics is not Bohmian mechanics. Note this 
Hamiltonian is unreasonable if X is position as this 
says we have a cubic dependence on momentum, which 
is not present in any natural Hamiltonians. However, 
if X corresponds to momentum and —Y to position, 
(which occurs in in Brown and Hilcy's [25] extension of 
Bohmian mechanics to include the momentum represen- 
tation), then this Hamiltonian is valid. 

As before if we discrctizc the value of X, 
the probability current again will be given by Eq. (3.14). 
However, the discretized version of the Hamiltonian in 
this case is 



(en|7J un i|em) = (iS n m+3 - iS n m - 3 - 3iS n m+i 

) + V{en)5 n , m , (3.28) 



(3.24) which results in 
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J nm (^) 



rim 



+3iS 

71,771—1 



(3.29) 



Taylor expanding this gives a rather large expression, but 
since (ip n (t)\i[>n(t)) is a real number and 1/e 3 3> 1/e 2 in 
the e -t limit we can ignore all orders of the Taylor 
expansion. This allows us to write 



J nra (^) 



e 3 

+35 n ,m-l]- 



^n,m+3 ^n,m— 3 3(5 n , m -|_i 

(3.30) 



Since (V'nMIV'raW) is always positive, for k > the tran- 
sitions rates defined in Eqs. (2.5) - (2.8) become 



T( n -3)n(t) = 



Ti 



(n+3): 



,(*) 



T( n -l) n {t) 
,(*) 



r, 



(n+l)' 



K 

3k 
0. 



(3.31) 
(3.32) 

(3.33) 
(3.34) 



These transition rates imply that in an interval dt the 
average displacement dx will be 

M[dx] = 3eT( n+3)n dt - eT^^dt = (3.35) 

and the dispersion will be 

M[dx 2 } = 9e 2 T (n+1) „ + e 2 dtT [n _ 1)n dt = Undt/e (3.36) 

which diverges as e — > 0. Thus the continuum limit does 
not exists. This implies that Brown and Hiley's [25] ex- 
tension of Bohmian mechanics to include the momentum 
representation (the momentum projector is the preferred 
projective measure) is not the continuum limit of Bell's 
modal dynamics. It is possible that a different choice for 
Jnm{t) (and T nm (t)) would allow their equations to be 
derived, but that is beyond the scope of this paper. 



B. The Velocity Operator Technique 

In the above section we have demonstrated that when 
using the Bell solution for the transition rates, modal dy- 
namics for some continuous properties only reduces to a 
deterministic theory (apart from a random initial con- 
ditions) if the Hamiltonian is at most quadratic in the 
conjugate variable to the property. If the property is po- 
sition then this deterministic limit is Bohmian mechan- 
ics and the trajectories are then found using Eq. (3.5), 
which requires calculation of the probability current den- 
sity, Jk({xj},t). Here we present an alternative to this, a 
method to calculate Vk({xj},t) directly. We assume that 

„ - ; f) Mm)\{xj})({xi}\Mt)m))] ( o o 7) 

- mt)\{*i})({*j}Mt)) . ' 1 J 

where Vk(t) is the k component of the velocity operator. 
This operator is defined as 



[Xk, H un i(t)). 



(3.38) 



To show that this docs give the same trajectories 
as Bohmian mechanics, we note that the numerator of 
Eq. (3.37) should be Jk({x 3 },t) by definition. Now using 
the Hamiltonian defined in Eq. (3.6), the velocity opera- 
tor is 



t> fc (t) = B k (t) + B\{t) + 2C k {t)Y k + 2Y k Cl(t). (3.39) 



This results in the following velocity field 



V k ({Xj},t) 



m{x 3 },t)\i;({x 3 },t)) 



-Re 



$({ Xj },t)\{2MBk({xj},t)} - 4iRe[C k ({ Xj },t)]d Xk 



-2id Xk [Ci{{x 3 },t)]\mx 3 },t)) 



(3.40) 



Comparing this with Eq. (3.9) we see that the numera- 
tor is indeed Jk{{x 3 },t). This completes our proof that 
our velocity method does generate the same trajectories 
as Bohmian mechanics for Hamiltonians of the form dis- 



played in Eq. (3.6). However, by extending this argument 
to higher orders it can shown that our velocity method 
does not agree with Bohmian mechanics. That is, this is 
another example of a method that only works for Hamil- 
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tonians that do not contain terms of order Y£ and higher. 



IV. SIMPLE EXAMPLE - HARMONIC 
OSCILLATOR 

A. Husimi POM 



where \(x + + x') / \f2) UYl \ is an a;-state (an eigenstate of 

x = (a + o+)/v^). 

Because x + and y~ can be well-defined simultaneously, 
we interpreted them as being suitable modal properties 
to represent simultaneously the position and momentum 
of the harmonic oscillator. This can be justified on the 
grounds that for \<j>) = |0) (a vacuum state) 



To illustrate modal dynamics for POMs, we investigate 
a simple model; a universe consisting of a one dimensional 
harmonic oscillator of frequency u. That is, 



H un i(t) 



tdcra, 



(4.1) 



where at and a are the creation and annihilation opera- 
tors of the harmonic oscillator respectively. 

The preferred POM we consider is that of Ref. [22], 
which has POM elements (or effects) given by 



F a = — \a){a\d 2 a, 

7T 



(4.2) 



where a\a) = a\a). The preferred observable (one of the 
many) we associate with this POM is 



A = {{a,-\a){a\d 2 a)} 



(4.3) 



The continuous value a = x + + iy~ is a complex num- 
ber representing a point is phase space (x + ,y~). In the 
orthodox theory this POM corresponds to a measure of 
both position and momentum with minimal additional 
uncertainty [18]. 

Before analyzing the modal dynamics that corresponds 
to this POM, we need to define a few operators that act 
in the enlarged Hilbert space K.. (The auxiliary Hilbcrt 
space Haux is assumed to be a single harmonic oscillator). 
We define 



- = [a + a* + b + P]/2, 
= [6 + fit - b- 6 f ]/2, 



= [—ia + id) - ib + ib ] ]/2, 
= [-ia + ia) + ib-iP]/2, 



(4.4) 
(4.5) 
(4.6) 
(4.7) 



where b and w arc annihilation and creation operators 
which act in H aux . These four operators obey the com- 
mutator relations 



[x + ,y + ] 
[£ + ,y~] 



F ,y 



+i 



0, 



(4.8) 
(4.9) 



thus x + and y have the same eigenstates, which we 
denote as \x + ,y~). They are given by 



r dx' 


x + — x' \ 


1 


\[2 ' aux 



V2 



iy x 



(4.10) 



mt)\y-y~mt)) = 



(*(t)|*|*(f))/V2, (4.11) 

(*(t)\Y\*(t))/V2, (4.12) 

(*(t)|l 2 |*(t))/2 + l/4, (4.13) 

(*(t)|F 2 |*(t))/2 + 1/4, (4.14) 



where = | x I , (t))|0) is the guiding wave for K, and 

\^(t)) is the solution of the Schrodingcr equation. Thus 
these operators have essentially the same statistics as the 



position (X) and momentum (Y = (- 



)/V2) 



op- 
erators in the classical limit (as the 1/4 term becomes 
negligible) . 

In K, we can rewrite H un i as 

-Huni^laux = Uj[x + X + + X~ X~ + j) + y + + jj~ jj~ 

+2x + x- + 2y+y- - 2]/4. (4.15) 

To define the modal dynamics in K, we use Naimark 
theorem and a Naimark projector \<p) = |0) to extend 
the POM elements defined by Eq. (4.2) to the projector 
y~) y~ \dx + dy~. That is 



-\a)(a\d 2 a = (0\x + , y~) (x+ , y-\0)dx + dy- (4.16) 

7T 

(see appendix A) . The observable A becomes a property, 
which in matrix notation is given by 



A = 



(x++iy )\x+,y )(x+,y \dx+dy . (4.17) 



Since the preferred projective measure in /C, 
{Tl(x + ,y~)dx + dy~ = \x + ,y~)(x + ,y~\dx + dy~} forms 
a complete orthogonal set and Eq. (4.15) contains no 
cubic or higher order terms involving x~ or y + , the 
results of Sec. Ill B are applicable to this paper. That is 
a deterministic differential equation for the values x + (t) 
and y~ (t) can be determined. 

Using Eq. (3.38) with the Hamiltonian Eq. (4.15) gives 
the following two velocity operators 



6_(t) 



ui 



2 2 

Substituting these into Eq. (3.37) (with |^(f)) 
gives 



(4.18) 
(4.19) 

!*(*))) 
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v+[x + ,y ,t) 
v_(x + ,y~,t) 



u;Bc[-im)\x + , y~)d x+ (x+,y-m))] , u 
2{*(t)\x+,y-}(x+,v-\$(t)) 

LjRc[i($(t)\x + ,y-)d y - (x+,y-\$(t))] 

2($(t)\x+,y-)(x+,y-mt)) 2 



Thus the differential equations are 
d t x + {t) 

dtv~{t) 



+ _ u;R C {-i(<S>(t)\x+, y~)d x+ (x+ , 



2($(i)|x+,y-)(x+,y-|$(t)) 
wRe[i($(t)|ar+ jtf -)a s - »-)$(*))] 



2($(t)|a;+,y-)(a;+ )V -|$(t)) 



x+=x+(t),y-=y-(t) 2 
x+=x+(t) l2 /-=s/-(t) 2 



, At — ^ 



Using the fact that (see appendix A) 

(x+,y-\$(t)) -- 
the partial derivatives can be written as 

d x+ (x+,y-\$(t)) = ex P[-^ +2 +^' 2 )/ 2 ] y m {x+ Zjffrl {mm)) - x +(x+,y 



>.! 



d v -(x+,y-\$(t)) 



(4.20) 
(4.21) 

(4.22) 
(4.23) 

(4.24) 
(4.25) 



-t — > m == (m wm) —y vxr.y \<P(t)). (4.26) 

. ITT ^ 



which allows us to write 



d t x+(t) = +^y-(t) + ^lm[ Xil (x + (t),y-(t))}, (4.27) 



dty~(t) 
where 



-|s + (t)-|Re[^(s + (f),ir(t))]> (4-28) 



x*(x + (t),ir(t)) 



T, m m(x+ - iy-) m - l {m\y{t))/yfm\ 



E m O+ ~ %-) m (m|*(t)) /Vm! 

(4.29) 

Thus the differential equation for a(f) is 



<M*) = -y «(*) ~ yXv(* + (*)> iT(*)) ( 4 -30) 



|n). Substituting this into Eq. (4.29), 



gives 

X*(:c + (t), »"(*)) 



n— 1 



n[a: + (t)-iy-(t)] 
[a;+(t)- l2 ;-(t)]™ 
na(t) 



|a(i)| a 

Thus 

*«(*) = -y(i + |; 

This has the solution 

a(t) = a(0)e-*"'* 



(*)l s 



a*(i) 
(4.32) 



a(i). (4.33) 



(4.34) 



where oj' = lu\1 + n/\a(0)\ 2 J /2. This solution and all 
subsequent solutions are discussed in Sec. IV C. 



1. When \^(t)) is a number state 

Lets first of all consider a number state \n), as the 
initial condition for |\£(0)). Then by the Schrodingcr 
equation 



dt|*(t)> = -«wota|*(t)), 



(4.31) 



2. When \^(t)) is a coherent state 

If we assume that initially the system in is a coherent 
state |*(0)) = \(3), then by Eq. (4.31), 

an p — iwnt 

|vl/(t))=cxp(|/3| 2 /2)y ; r= \n). (4.35) 
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Substituting this into Eq. (4.29), gives 



1. When is a number state 



,J2 m m ( a *l 3e ' iWt ) m ~ 1 /' m ^ 



E m (a*/5r Mt ) ra /m! 
= (3e~ iuJt . (4.36) 

Thus 

dta(t) = ~a(t) - y/3e-^*. (4.37) 

This has the solution 

a{t) = MO) - (3}e-' 1 ^ 2 + f3e~ luJt . (4.38) 

B. Bohmian Mechanics, The Position Projector 

In this section we consider the modal dynamics 
which describe a decomposition into position cigcnstates. 
That is the preferred projective measure is {n(x)dx = 
\x)(x\dx}. Since this is already a projector, there is no 
need to enlarge the universe, and as shown above (and in 
Ref. [4] and Rcf. [5]) the modal dynamics for this case 
is just Bohmian mechanics. 

In terms of X and its conjugate operator Y, the Hamil- 
tonian for the universe, Eq. (4.1), becomes H un i{t) 



H uni (t) = co(X 2 + Y 2 -l)/2. 



(4.39) 



Using our velocity operator technique it can easily be 
shown that the velocity field is 



v(x, i) 



LuRc[-i(y(t)\x)d x (x\y(t))} 
(*{t)\x){xm)) 



as v — loY . Using Eq. (3.5) this gives 
wIm[(*(t)|x)S x (s|*(*)>] 



d t x(t) = 



(*(t)\x)(x\*{t)) 



x—x(t) 



(4.40) 



(4.41) 



Since Y does not commute with X, we can not give 
both X and Y definite status (to give Y property status 
we would have to chose a momentum projective mea- 
sure as the preferred projective measure). However, as 
in Bohmian mechanics, we can define a momentum field, 
V(x,t), by 



y(x,t) 



Re{(*(t)\x)(x\m(t))} 

(*(t)\x)(xMt)) 
hn[(*(t)\x)d x (x\*(t))] 



(*(t)\x)(x\9(t)) ' ^ 
One interprets this momentum field as follows: if the 
system has the position x(t) then its momentum is 
y{x)\x=x(t)- With this and position x(t) we can define 
a phase point (x(t), y(x)\ x=x ^), which in complex nota- 
tion is written as 



a(t) 



x{t) + iy{x)\ x=x{t) 
V2 



(4.43) 



The factor l/y/2 is to scale this to agree with the preced- 
ing section. 



As before when we assume a number state initial con- 
dition the guiding wave at time t is = e~ tujnt \n). 
Substituting this into Eq. (4.41) and using 



J-ex P (-x 2 /2)^ 



H n (x) 



where H n (x) is a n th order Hcrmitc polynomial, 



(4.44) 



d t x(t) 



LuIm{d x [exp{-x 2 /2)H n (x)]/v / 2 n n\} 



cxp{-x 2 /2)H n {x)/V2^. 



= 0. 



x—x{i) 

(4.45) 



Using a similar argument and Eq. (4.42) we get 
y(x, t) x — x u\ = 0. That is, once the initial x(0) is picked 
from the quantum mechanical distribution is stays there 
for all time. In terms of the complex notation a(t) we 
get 



y(t) = x(0)/V2. 



2. When |\£(t)) is a coherent state 



(4.46) 



If we assume that initially the system in is a coher- 
ent state |*(0)) = then Eq. (4.35) in the position 
representation is 

!*(*)> = 6XP y 2) J dx' eMV2/3e~^x 

-P 2 e- l2ult /2) exp(-a;' 2 /2)|x'). (4.47) 
Substituting this into Eq. (4.41) gives 

u}lm{d x [cxp{-x' 2 /2) cxp(V2Pe- luJt x)}} 



d t x(t) 



exp 



(~x' 2 /2) exp(v / 2/3e- lwt x) 



)V2Im[/3e 



and using Eq. (4.42) we get y(x, t) x=x ( t ) = V / 2Im[/3 | 
Taking the derivative of this gives 



<ky(x,t)x=x(t) = -u>V2Re[(3e~ 



x—x(t) 

(4.48) 

, — iujt j 

(4.49) 



Thus 



d t a(t) = Lu{Im[(3e- tuJt ] - iRe[f]e- tuJt }} = -iv /3e- iu,t , 

(4.50) 

This has the solution 



t(t) = f3e- tult + (a(0) - (3). 



Classical Limit 



(4.51) 



The classical harmonic oscillator has the well known 
solution, this being 



a(t) = a(0)e" 



(4.52) 
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In Sec. IV A (modal dynamics when the preferred POM 
is the Husimi POM) we saw that in the enlarged Hilbcrt 
space, wc can define a(t) = x + (t) +iy~ (t). Now consider 
the classical limit. First consider the case of the number 
state Eq. (4.34). From the probability formula (2.22) it 
can be shown that |a(t)| 2 m n with high probability. This 
means that as n — > oo, a/ — > lj 1 reproducing the classical 
dynamics. When considering the case with an initial co- 
herent state we can similarly argue that (3 « o;(0) with 
high probability. Then in the limit \(3\ — ► oo, the dif- 
ference between the classical formula and Eq. (4.38) is 
negligible. 

By contrast, in the position case (Bohmian mechan- 
ics), for the first case (number state), a(t) = x(0) for 
all n, the dynamics is completely non-classical. How- 
ever, it can be argued that when wc consider the second 
case with an initial coherent state, again the difference 
between the classical formula and Eq. (4.51) is negligi- 
ble. The fact that in both Bohmian mechanics and the 
modal interpretation of Sec. IV A have a good classical 
correspondence for coherent state is not surprising, as the 
coherent state is a classical-like state. What is surprising 
is that it is possible to obtain classical modal dynamics 
even for a non-classical state, by using POMs. 

V. DISCUSSION 

In this paper we have extended the modal dynamics of 
Refs. [3, 4, 5, 7, 8, 9] to include the possibility of hav- 
ing a preferred POM. To do this we enlarged the Hilbcrt 
space from H un i to /C = H un i <S> H aU x- Once in this en- 
larged Hilbcrt space we used Naimark's theorem [18, 24] 
to define a preferred projective measure which is equiv- 
alent to the preferred POM. That is the statistics of an 
observable describe by Z = {(z n , F n (t))}, (where F n (t) 
is a POM element) is equivalent to the statistic of the 
property Z = {(z n , tl n {t))}- Here II(i) is a projector in 
K, and defines the property state | <!>„(£)). This state rep- 
resents the actual state of the universe and determines 
the present value, z n , of the property Z from the set of 
possible values. The evolution (jumping between z n ) is 
determined by the stochastic evolution of | <&„(£)), which 
in-turn depends on the guiding wave = \<f))\^{t)). 

\^(t)) is the standard quantum state found from the 
Schrodingcr equation and |</>) is a state defined in the 
auxiliary Hilbcrt space. 

To illustrate modal dynamics for POMs we considered 
a simple example: a universe consisting of a single Har- 
monic oscillator. To illustrate our new dynamics, we 
looked at the Husimi POM and compared the dynam- 
ics obtained to that which is obtained with the position 



projective measure (Bohmian mechanics). Since the first 
case corresponds to a POM we have to use our enlarged 
Hilbert space modal dynamics to develop the stochas- 
tic evolution equation for the value of the property (or 
equivalently the property state). For the Husimi POM 
we denoted this value by a(t). We find that by choosing 
a Naimark projector \<p) = |0), this property defines a 
point in phase space (a(t) — x + +iy~) that is effectively 
the position and momentum of the system. Investigat- 
ing two different initial conditions for |^>(0)), namely a 
number state and coherent state, we find that the differ- 
ential equation for a(t) for both cases has a classical limit 
which agrees with classical theories. When comparing to 
the modal dynamics for the position projective measure 
we find that when the initial state is the number state 
the dynamics are highly non-classical. Only for an initial 
coherent state (a classical-like state) can a classical limit 
can be obtained. 

In conclusion by extending modal dynamics to include 
POMs allows the ability to include overcomplete decom- 
position, like the Husimi POM [22, 23] and information- 
ally complete POMs [20, 21]. This may provide an an- 
swer to questions involving the quantum-classical limit. 
An interesting question for future work is whether the ex- 
tension of the Hilbert space is only a mathematical tool 
or whether there is some physical significance behind en- 
larging the Hilbert space. Apart from this fundamental 
question we intend to use this theory to explain diffusive 
non-Markovian stochastic schrodingcr equations (SSEs) 
[26]. In a recent paper we have shown that under the 
orthodox interpretation non-Markovian SSEs represent 
nothing more then a stencil which determines the state 
of a system at a particular time t, given that a measure- 
ment on the environment at that time yielded result z 
[27]. That is, non-Markovian SSEs can not be interpreted 
as evolution equations for the state of the system condi- 
tioned on the outcomes of some continuous measurement 
of the environment. We believe that under the modal 
theory it can be shown that non-Markovian SSEs are 
evolution equations for the system part of the property 
state when the environment property Z is given definite 
status. In [28] we have shown that when the preferred 
measure is position, non-Markovian SSEs do have this 
interpretation. To interpret other non-Markovian SSEs 
[26, 27], it is necessary to consider non-orthogonal de- 
compositions. 



APPENDIX A: PROOF OF EQ. (4.16). 
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To show that |a)(a|/7r = (0\x + ,y )(x + ,y |0) we consider just (0\x + ,y ) = |a)/y / 7r. Using Eq. (4.10), and the 
standard definition of a x-state Eq. (4.44), (0\x + ,y~) becomes 



(0\x + ,y-) 



dx' 1 



cxp 



/2tt 7rV4 

Defining X = (x + + x')/y/2 allows us to rewrite this as 



-(x+ -x'f 



dX 1 



■ cxp 



-{2x+-^2Xf 



/tFttI/4 

which with definition Eq. (4.44) can be expanded to 



V2 



\x) uni V2e*y-^ 2X - 



/dX 
exp(-a; + a; + + V2Xx + - X 2 + iy~V2X 



ly x 



H n (X) 



exp(-|a| 2 /2)E J ^exp(-X 2 )exp(2aX/V2-a 2 /2)^^l 



where a = x + + iy . Then using J2 m t m H m (x)/m\ = exp(2ta - t 2 ) this can be written 



<\x+,y-) 



exp(-|a| 2 /2) 



7T 



E 



/ 2™m!V2"?i! 



dXeM-X 2 )H rn (X)H n (X) 



and since J dX exp(— X 2 )H m {X)H n {X) = 2 n nly/nS nm this becomes 



(0\x + ,y-) = 



cx P (-|a| 2 /2) 



Thus \a){a\/ir = (0\x + , y~) (x+ , y~\0). 
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